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Summary~— The movement inside rising bubbles has been investigated by visual and indirect 
means. Gas bubbles of smaller diameter than 0-03 em in water behave as rigid spheres, and 
above this diameter a vigorous toroidal circulation sets in, which persists throughout the rise. 
The transition from rigid to circulating conditions has been determined for several liquids, 
and it has been found that the critical radius of transition is not in accordance with the formula 
put forward by Boxp and Newton. 

The velocities of bubbles in water have been measured and the discrepancies between 
the data in the literature have been considered in terms of four variables : temperature, rate 
of formation, wall effect and contamination from surface active agents. 


Résumé—Les auteurs étudient par lobservation visuelle ou par des méthodes indirectes le 
mouvement a lintérieur des bulles gazeuses au cours de leur ascension. Pour des diamétres 
inférieurs a 0,03 em., les bulles gazeuses se comportent dans Teau comme des sphéres rigides. 
Pour des diamétres plus grands, il s’établit une active circulation toroidale qui persiste tout 
au cours de lascension, 

Les auteurs ont déterminé pour différents liquides le seuil de passage de l'état rigide a l'état 
circulatoire : ils ont trouvé que le diamétre critique de transition ne correspond pas aux valeurs 
déduites d'une formule proposée par Boxp et Newron. Ils ont mesuré les vitesses des bulles 
dans l'eau et les écarts que présentent les différentes valeurs trouvées dans la littérature ont été 
examinés par rapport aux quatre facteurs : température, vitesse de formation, effet de paroi, 
contamination de la surface par des agents tensio-actifs. 


INTRODUCTION 

The velocities of rise of a range of sizes of gas 
bubbles in water were required for mass transfer 
calculations, and the internal motion in bubbles 
was investigated in order to interpret the result¢. 
A recent survey [8] of the available data on 
velocity of rise of bubbles revealed serious 
disgreements in the diameter range from 0.06 
to 0-1 cm. 

In fluid droplets or gas bubbles the existence 
of a regular current of fluid up the central 
axis and down the sides of the drop had been 
predicted [15], [29] and demonstrated [11], 
[12], [13] to be due to the viscous drag of the 
outer fluid. The velocity of rise with internal 
circulation as postulated by Hadamard [15] 


should be 1-5 times that of bubbles in which no 
circulation occurs~that is of bubbles which 
behave as though they had the rigid surface of 
a solid sphere. This ratio of the actual velocity 
to the theoretical velocity calculated from 
Stokes’ law can be used as an indication [4] 
of the existence of circulation, and is denoted 


as ; in the equation 


V = a) 


where Vis the velocity of the bubble; Ap the 
density difference between phases ; g the accelera- 
tion due to gravity and ,« the absolute viscosity 


of the outer phase. The value of i of 1-5 for a 
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circulating bubble and 1-0 for a non-circulating 
bubble holds only in the range of ReyNo.ps’ 
number where Stokes’ law is valid. SroKkes’ 
law for rigid spheres or for gas bubbles or liquid 
droplets behaving as rigid spheres is normally 
only valid for ReyNoips’ numbers appreciably 
below 1. If the velocity of circulation in gas 
bubbles or liquid droplets is as postulated by 
Hapamarp, that is equivalent to the rate of 
fall or rise of the sphere, then the stream lines will 
be similar to those in the Stokesian region which 
therefore should extend considerably higher than 


1 
1. Thus we can use i at much higher ReyNoups 


numbers than those in which it is normally 
assumed that Stokes’ law holds. It must also 
be assumed that there is no distortion of the 
droplet or bubble shape from a sphere although 
in fact such distortion does occur at relatively 
low ReyNo.ps’ numbers for air bubbles. 

Bonp and Newton [4] determined the transi- 
tion diameter in four systems : air in waterglass, 
air in Golden Syrup, (sugar solution), mercury 
in Golden Syrup and water in castor oil. An 
attempt was made to characterise the transition 
region by the dimensionless group RT W which 
reduced to a critical radius defined as : 


r = (T,Apg)! 


OBSERVATION OF CIRCULATION 


(2) 


Motion inside gas bubbles was investigated by the 
addition of ammonium chloride fog during for- 
mation. The most suitable method was to bring 
together streams of air containing ammonia and 
hydrogen chloride with one component in excess. 
Jets of different shapes, both upright and inverted, 
were used and demonstrated that circulation 
was not due to the energy of the gas leaving 
the jet. 

The liquids studied were glycerol, white 
hydrocarbon oil and water. In liquids of viscosity 
less than 1 poise it was not possible to observe 
movement in very small bubbles (smaller than 
0-35 cm dia.) or in any very rapid bubbles. 


Glycerol 
Circulation was observed at diameters between 
0-44em and 1-2cm. The velocities of bubbles 
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‘able 1. Divergence from law of bubbles in 
glycerol. 
Bubble Bubble 
dia. velocity Wall Lk 
em, em see Effect (corrected) 
Air with Fog 
OOS 105 139 
O-41s 0-046 13% 1-22 
0-445 0-057 16° 
O67 2-09 23° 
12 37 0-68 45°, 
Carbon 
445 mow 7 % 
O-359 874 O-026 
0-363 ONTS | 0026 9°, 
0-478 14500 0-057 10% 1-53 
0-488 0-061 11%, 1-33 
0-610 224 14%, 
0-640 251 0-183 13%, 1-58 
liv 
0-235 O3TL 0-007 | 50% 
0-403 52% 138 
O-364 0-93 8-7 % 1-33 
109 0-048 92% 
| 10-6 % 1-50 
0-458 163 0082 | 
0-500 1-85 | 20% | 1-46 
0-530 203 | 123% 144 
O580 | 262 O165 144% 155 
0666 320 65% 148 
0-732 350 O280 180% 
0-760 38 otis «200% 
1-03 6-67 O72 | 248% 
7-78 os | 290% 
1280 | 310% | 122 
1:33 366% | 108 
Sulphur diovide 
622 68%, 144 
0-299 0-767 0-023 72% 1-61 
0399 0-051 95% | 1-36 
O46 0-054 100% 
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+—WATER 
10 
REYNOLDS NO, 
Fig. 1. Divergence of the rate of rise of bubbles from Strokes’ law, 


with and without fog were measured, and 
expressed in terms of i (Table 1). Fig. 1 shows 
the results in dimensionless form. At diameters 
below 0-20 cm (Ny, = 0-002), the value of ; is 1-0 


which indicates that the surface is rigid ; above 
0-30 em (Np, = 0-023) there is agreement with 
Hapamarp’s correction, thus confirming the 
visual observations. The transition diameter 
between the two regimes was defined by Bonp 


and Newron [4] as the value of d when 7 = 1-25 


and for this system it is 0-25em + 20%, 
The full curve in Fig. 1 through the points 


for glycerol falls sharply from | = 15 at Np, = 


0-2 and becomes identical with the curve for 
rigid bubbles, the lower curve, at N,, = 6-0. 
Thus the value of i cannot be used to detect 


circulation in bubbles with Reyno.ps’ numbers 
of 6 or over, and, as mentioned later, the bubbles 
are no longer spherical. 


White Oil 
Streams of bubbles containing fog were injected 


into a hydrocarbon oil (viscosity 1-98 stokes 
at 20°C) at a frequency of 0-1 to 3 per second. 
Circulation was observed over the diameter 
range 0-3 to (Np, = 0-5 — 1-8). The 


values of i were all greater than 1-5 owing to 


the rapid rate of formation. Values as high as 
2-5 were obtained at frequencies higher than 3 
bubbles per sec. 

A further series of velocity determinations 
was made with single bubbles with diameters 
from 0-04cem to 0-25em. The bubbles were 
formed without fog, and were measured in a 
thermometer capillary. The results given in 
Table 2 agree with Hapamarn’s equation at 
all diameters. The transition diameter is there- 
fore less than 0-04 cm. 


Water 
The transition diameter for water could not be 
observed directly, but could be calculated 


ond from date 


on mass transfer from bubbles. Below 0-015 em 
dia. (Ng, < 1-1) the value of i is 1-0 (within 8%) 
and therefore the bubbles rise as rigid spheres. 


approximately from values [1] of 
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Table 3. Transition radius, em. 


Table 2. Hydrocarbon white oil. 
| | | | | | 
Bubble — Bubble Bond 
dia. velocity (Wall Wk System Observed ry Reference 
em. | em/see | Np, Effect (corrected) radius | radius | 
0-036 0-063 O-00124 | 5% 1-75 Air in waterglass Onl 025 Bonn and 
0-064 O-137 0-00475 10%, 1-25 Air in syrup O11 0-25 
0-066 O-195 0-00695 1-62 Mercury in syrup O20 [4] 
0-085 0268 13%, Water in castor oil O77 0-70 
O-116 O476 19%, 
O515 0-827 20%, 132 Air in white oil 
0-125 | 0-592 Oow 21% 155 Air in glycerol O125 | 0225 | Present 
O-254 2-08 28S 43% 1-57 Air in water 0020 O27 Work 
O-254 1-88 257 43% 142 Air in water + 2% 
Cellofas 0-25 0-25 


The rate of mass transfer from bubbles of 0.03 em 
dia. is reported [26] to be 45cm, hr, and for a 


circulating bubble is hr. Circulation 
therefore starts between 0-015 and 0.03 cm dia., 
but it does not reach the maximum speed until 
0-25 em dia. [16). 

Circulation in bubbles can be observed if the 
viscosity of water is increased to 115 c.p. by 
addition of 2%, of Cellofas B, a hydroxylated 
cellulose. Circulation was observed with bubbles 
larger than 03cm dia., and the transition 
diameter was between 0-3 and 0-5 em as shown 
both by mass transfer data and by values of 


i: A similar value for transition diameter was 


observed for a 1% solution having a viscosity 
of 30 ¢.p. 

It is shown later that circulation at all dia- 
meters up to at least 0-6cm can be stopped by 
very small quantities of surface active materials 
adsorbed on the bubble. It is inferred that the 
“ condensed solid" type of film [17] is formed, 
making the surface rigid. 

Data from the literature 

ARNOLD [2] reported the velocity of bubbles in 
olive oil and linseed oil, from which values of 7 
were calculated. His results show that in olive 
oil, circulation starts at 0-Olem dia. and is 
fully established at 0-07 em dia. Above this dia- 


meter there is agreement with the Hadamard 
equation (Fig. 1). Bubbles in linseed oil have a 


| ALLEN [1] 


Air in aniline > 0-05 O42 
Air in olive oil O04 | 0-173 ARNOLD 
[2] 
Air in linseed oil > O15 — ARNOLD 
{2} 


value of 1-0 (within 3°) until inertia effects 
become important at Np, = 0-2. 


Discussion 
In Table 3 are given the observed transition 
radius and that calculated by and Newton's 
equation (2). There is a marked divergence, 
between these values and only for water in 
castor oil and for air in a Cellofas solution is 
there good agreement. 

The failure of equation (2) to predict the 
transition diameter implies that other variables 
exert an influence on circulation, and the move- 
ment of the surface must be considered in more 
detail. The surface of the bubble moves outwards 
and downwards from the forward stagnation 
point, thus forming new surface. The energy 
required per second to form the surface is 
2ar Tu ergs, where u is the tangential velocity 
of the surface at the ‘‘ equator” of the sphere 
and is the interfacial tension, dynes /cm. 
This energy is given up from the rear of the bubble 
by the closing up of the surface, and there may 
be a net loss of energy if the interfacial tension 
in compression is less than in stretching. Bryn 
[5] refers to this difference as the relaxation 
effect. 
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Table 4. Table 5. Rate of rise of bubbles contaminated with 


Diameter | Velocity Temperature | 
em/sec. | Cc | Nae Diameter Velocity Temperature) Time since 
——_—_| — —- em | em/sec c previous 


Air bubbles in Voom dia. column bubble 


* Vaseline.” 


cm 


0-0186 145 14-6 Air bubbles 
0-0258 1-96 | 146 
0-0307 27 8 | 146 o-s16 | 30min. 
0-0313 282 ame. 
0-0821 282 13-5 10 min. 
0-0326 3-17 143 0-282 12 hr. 
0-0357 3-50 146 0-291 72 hr. 
0-0357 3-45 140 = 
0-0402 4-02 146 | | 20 min. 
0-0417 372 | 148 0-578 min. 
0-0428 34 146 0-423 | 5 min. 
0-0490 5-43 10 0-440 20 min. 
0-0512 5-75 140 | | 
0-0522 5-65 12-5 owe ain. 
0-455 , 30 min. 
Air bubbles in 2em dia. column 0-484 
| 0-525 
O-0126 O77 | 150 0-518 
OBS 0-534 
0-0205 148 | 0-550 
0-0252 2-26 
0-0336 2-65 0-575 
00-0359 3-22 0 0-619 
O-0572 a4 50 0-655 
0-678 
Air bubbles in dia. column O-745 
0-307 246 18-0 0-783 
O-314 246 18-0 
0-400 24-0 18-0 
0-425 23-4 160 0-253 
0-441 22-2 18-0 0-374 
0-460 21-6 18-0 0-435 
0-483 16-0 0-528 
0-492 23-2 16-0 0-634 
0-499 22-7 18-0 
0-504 2240 180 
0-675 16-0 


1560 0-573 


te Thus any change of surface tension with time 
would be an important factor in determining 
O16 | se 178 the start of circulation, since if the work done 
0-2 21-7 17” on the bubble by the skin friction is less than 
0-23 | 2-4 170 any nett loss of surface energy, the circulation 
0-25 22-6 170 5 cannot persist. 
om ae ra The transition diameter can only be completely 
in ae ai correlated for all systems if terms for the intensity 
0-32 get IT” of skin friction and for the change of interfacial 
tension with time are included in the expression. 
0-453 22-6 17-0 More data are required on transition diameters 
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for pure systems and for those containing sur- 
face active agents. 


MEASUREMENT OF VELOCITY 

The velocity of small bubbles was measured 
with a stop watch over a rise of 73 em. Columns 
of 10cm and 2m dia. were used, containing 
water saturated with air. Bubbles of diameters 
down to 0-0125 em were formed singly in a ther- 
mometer capillary [7] of 0-005em dia. Below 
0-010 cm dia. the errors of measurement are 
greater than 10°, by this method owing to the 
curvature of the meniscus of the confining 
liquid. 


OSCILLATION AND CHANGE OF SHAPE 


Bubbles smaller than 0O-lem diameter are 
approximately spherical and rise in straight 
vertical lines. At a critical diameter reported 
by Horrer [18] as 0-10em and by Bryn [5] 
and Lucusincer [20] as 0-15em, the path of 
rise becomes a helix. The pitch and the amplitude 
of the helix are both approximately two diameters, 
The bubble is flattened into a lenticular shape, 
and moves in the direction of the major axes. 

The dimensions of the helix increase with the 
diameter up to 0-5em diameter, when periods 
of straight vertical rise interrupt the helix. 
Bubbles bigger than 0-8 cm diam. nearly always 
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DRAG COEFFICIENT 


REYNOLDS NO.. Np, 


Fig. 2. 


Correlation of drag coeflicient and Reynolds’ 


No. for bubbles in water. 


Larger bubbles of 0-15 cm to 0-8 em diameter 
rise too rapidly for accurate timing by stop 
watch. A cinephotographic method was used, 
having a repeatability of 1°,. The cine camera 
was brought to focus through a system of mirrors 
on two marks 56cm apart on the column. The 
camera speed was found to be constant within 
1% over both long and short periods, at a speed 
of 18-4 frames sec. Scrupulous care was necessary 
in the cleanliness of the apparatus to obtain 
consistent results. 


rise in straight lines. Their shape is hemispherical 
in front and flattened at the rear. Changes of 
shape are confined to irregular swinging of the 
rear part. Bryn and Hoerer report similar 
observations, but Corrock and MEIKLEJOUN 
[7] recorded that when bubbles rise in series the 
times of rise are erratic. 

Some bubbles in our work did not rise in the 
usual helical manner, but followed a zig-zag 
path, moving bodily from side to side in one 
vertical plane. This type of oscillation was found 


10 100 1000 1O000 
| 
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to be due to contamination from Vaseline 
used as a lubricant on the plug of a wide-bore 
glass tap. The bubble was held below the tap 
before being released into the timing section 
of the column. During this residence time, 
surface active constituents of the Vaseline 
were adsorbed onto the surface of the bubble, 
and mass transfer data [16] later showed that 
the surface was being prevented from circulating 
by a hydrocarbon film. 


EXPeRIMENTAL Resutts 


Table 4 gives the measured velocities of single 
bubbles in water, and Fig. 2 compares them in 
dimensionless form with the drag curve for solid 
spheres [6]. Bubbles smaller than 0-1 em dia. 
(Np, = 100) have about the same drag coeflicient 
as solid spheres. At higher ReyNoups’ numbers, 
the drag is less than for a solid sphere, even 
though the distance to be travelled along the 
helical path is greater than the vertical distance. 
This is because the bubble is flattened, and 
presents a considerably smaller frontal area [22] 
in the direction of motion than a sphere of equal 
volume. 

The velocities of bubbles whose surfaces are 
contaminated with surface active molecules 
(Table 5) do not show these low drag coetlicients 
at Np, = 200 to Np, = 1,200. Fig. 2 shows 
that the drag is greater than on either solid 
spheres or clean bubbles. At very high Reyno.ips 
numbers the drag coefficient becomes constant 
[9] at a value of 2-8. This value is higher than 
for a spherical body of the same volume because 
a greater frontal area is presented in the direction 
of motion. This difference in the direction of 
clean and contaminated bubbles has been referred 
to previously. 

The graphical presentation of data as a logari- 
thmic plot is not altogether satisfactory, since 
the accuracy of values taken from the curve is 
poor. A linear plot of velocity against diameter 
has been preferred in this work. Fig. 3 shows 
the experimental data for clean bubbles, and 
Fig. 4 the effect of contaminants on the rate of 
rise. 

The smallest bubbles are in agreement with 
the results of ALLEN [1] and Corpock [7]. The 
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20F 
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BUBBLE DIAMETER, CM 
Fig. 3. Velocity of gas bubbles in water. 


velocities are represented by the equations 
V = 108d for 0-01 <d <004em (8) 
V = 108d for 0.04 <d<0-20cem (4) 


In the diameter range 0-2 em to 0-6 em figures 
in the literature vary widely, since it is in this 
region that the four major variables 


(a) time interval between bubbles 
(b) surface contamination 

(c) temperature 

(d) wall effect 


have caused wide discrepancies. The effects 
of these variables are discussed below. 


Time InrervaL Berween BuBBLes 


The interaction between successive bubbles 
rising in the same vertical path has been calculated 
[3] for the Sroxes’ law region. At high ReyNoups’ 
numbers, the data available (Fig. 5) show that 
when bubbles are separated by less than 3 seconds 
the velocity is appreciably higher than for the 
single bubbles shown in the lower curve. The 
difference was found to be as much as threefold 
[8] at frequencies above 1 per second. 

Miyact quotes results that differ markedly 
from the present data, and although he states 
that formation was “ very slow,” it is evident 
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Fig. 4. Velocity of single bubbles with surfaces contaminated by Vaseline. 
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Fig. 5. Velocity of streams of bubbles rising in water, compared with that of single bubbles. 
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from his description of the photographic recording 
that bubbles were in fact formed at approxi- 
mately 1 per second. 

Coprock and MeIkLesoun [7] formed bubbles 
at one per second, and the velocities were found 
to be from 2% to 30%, greater than those of 
single bubbles. Accurate detailed study showed 
that the velocity over 10 cm. intervals varied by 
as much as 20%, from the terminal velocity 
over 200cm. The variation was random, and 
was not observed with single bubbles. The 


the frequency was about 1 per sec. Other workers 
[8], [19] [23], [24] have demonstrated the higher 
velocity of streams of bubbles at diameters 
above 0-5 em. 

It can be concluded that the error due to in- 
teraction of successive bubbles is less than 1% 
only when the interval between them is greater 
than 3 seconds, 


SuRPACE CONTAMINATION 
The cause of the lower velocity of contaminated 


d 
be 

O 
° 
w 
> 
w 7 
SF 

Ol O72 o4 Os 


BUBBLE DIAMETER, CM 


Fig. 6. Comparison of the data of Coprock and 
MrikLesoun with the curve for contaminated bubbles 
(from Fig. 4). Circles are data of Coprock and MEIKLE- 
Joun and the curve and line are taken from Fig. 4. 


irregularity of motion in streams of bubbles may 
account for many of the discordant results in the 
literature. 

Guyer and Pristerx [14] formed bubbles of 
carbon dioxide in water in a 3 cm. dia, column 
at a frequency of 8 per sec. The velocity is from 
50%, to 12% greater than that of single bubbles. 

Bryn [5] and Lucusincer [20] do not report 
the precise rates of formation. In both cases, 
bubbles were formed at a glass jet from a con- 
tinuous supply of gas, and it must be assumed that 


bubbles as shown in Figure 4 lies mainly in the 
shape of the oscillation and the greatest difference 
in velocity lies in the region where oscillation 
is most pronounced. The length of a helical path 
is greater than that for a zig zag rise of the same 
total height. The actual point velocities of 
contaminated and clean bubbles must therefore 
be even more widely different than the overall 
velocities. This appears to be due to the fact 
that in helical rise the bubble never comes to 
rest, whereas in a zig zag it “stalls” twice in 
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every oscillation. Skin friction is negligible 
compared with form drag at ReyNo.ps’ numbers 
over 100, but cessation of circulation has an 
effect not only on skin friction but will also 
cause some increase in form drag. 

The data of Coprock and Meikiesoun for 
single bubbles formed slowly at a jet are com- 
pared with the present data in Fig. 6. It appears 
that some contamination was occurring in their 
work during the residence on the jet, and possibly 
lubricant from the controlling needle valve 
travelled along the metal surfaces to the jet. 
The scatter of the experimental data is far 
greater than the error of the method, and may 
be the result of different degrees of contamination 
due to different times of contact with the jet. 


TEMPERATURE 


The present experiments covered the range of 
temperatures from 14°C to 20°C, and no variation 
of velocity was observed for either clean or 
contaminated bubbles from 0-2 cm to 0-8 em. dia. 

Data presented by Persies and Garner [27], 
however, shows velocities at 0-21 em and 0-32 cm. 
dia. to be 40°, and 14%, higher, respectively, 
than the present results. The temperature is 
not reported, but from the values of viscosity 
and surface tension used it appears to be 33°C, 
The differences cannot be ascribed to error, 
since a stop-watch can be read to within 10°, in 
a measurement of 2 sec. The frequency of forma- 
tion was reported to be 0-1 to 0.06 bubbles per 
sec., and if this applies to the two smallest 
bubbles then the difference must be due to 
viscosity. At diameters above 0-4 cm the velocity 
does not vary with small changes in temperature 
and viscosity, and by applying the appropriate 
wall correction, the data of Peesies and GarBer 
are in agreement with the majority of workers 


(5), [7], [9], [22], [28). 
Wa. Errectr 

The motion of spheres in the region where 
Stokes’ law applies is reviewed extensively by 
Barr [3]. Corrections for the presence of con- 
fining walls have been derived by LapeNBERG 
and by Faxen [3]. The Lapennerc correction 
is given in equation (5). 
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Table 6. 


Effect of confining walls on the velocity of spheres 


and bubbles. 


| 


& 


1810 
0-70 2110 
O71 2720 


O17 


0-35 

0-135 
0-225 
0-25 
Oo 
0-358 
0-154 


0-49 


System and 
Author 
Steel spheres in 
glycerol 


LuNNOoN [21] 


in water 


Air in water 
Corrock and 
Meikiesoun [7] 


Air in water 
| 


and 
Garser [27] 

Air water 
LucHSINGER (20) 


Air water 


O'Brien and 
[23] 
Air water 
Oyama and 
Iwase [25] 
Steel in 
water 


LuNNon [21] 


ad | aD K | Np, 
026 O11 0-20 | 
O38 O17 295 
0-5 0-22 0-62 | 
3 | 066 1-37 | 
om » 
004 O62 4 
005 0-025 | 
006 0-08 
0-07 0-085 

009 0-045 105 | 3 
0-10 | 0-05 Present: work 1954 
Be 
0-706 0-27 | 
O32 

108 | 

115 O44 | 

20 066 

20 0-33 5000 

28 0-18 
— |} | —_ 

20 O-4 6000 | 
| 

_| 

O64 | 0-128 7600 

0-96 | 0-185 13200 

128 | 25000 

1-60 0-308 25600 | 

1-98 | | 1400 

O64 7400 

0-96 0-326 12300 | 

128 | 37000 | 

160 0542 040 | 21000 

198 06651 24800 | 

10 
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V =v(14+21d/D) (5) 
where 


V is the velocity in an infinite fluid 

v is the velocity in a cylindrical container 
D is the diameter of the container 

d is the bubble diameter. 


The equation applies for values of d D below 
0-1 and for values of Ny, below 0-2. At high 
ReYNOLDs’ numbers there is satisfactory 
correlation. The data of several workers are 
summarised in Table 6 where the LapENBERG 
equation has been used as a basis of comparison 
because of its simplicity. The coefficient of the 
term d/D is designated K. 

The conclusions from Table 6 are that as 
Np, increases, K decreases: but that for any 
given column diameter A increases with Np,. 
The error introduced into determinations of 
velocity by the wall effect is less than 1°, if the 
column is more than 100 times the diameter of 
the bubble. 


CONCLUSIONS 


1. The velocity of a single bubble in water 
has been determined with an accuracy of 
2%. 

2. The four major sources of divergence from 
the true velocity and the limits of their 
effectiveness have been investigated. 

8. The existence of regular circulatory currents 
in gas bubbles rising through three liquids 
has been demonstrated. 

4. The diameter at which circulation 
commences in different liquids has been 
listed, and the Bond criterion has been 
shown to be unreliable for many of these 
systems. 


NOMENCLATURE 
d = diameter of bubble, em 
= gravitational constant, cm /sec®. 
K = constant in the Ladenberg equation (5) 
Np, = number, 
r = radius of bubble or drop, em. 
T = interfacial surface tension, dynes /em. 
u = tangential velocity of bubble surface, cm/sec. 
v = velocity of rise of bubble, cm/sec. 


Ay = density difference between phases, g/ce. 
» = viscosity of continuous phase, g/cm. sec. 


bubble velocity divided by Stokes’ law velocity 
p 
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Summary— Experiments and a study of the literature confirm the general equation u = Ve? 
where wu is liquid velocity above the bed, V' the free falling velocity of the particles in water and 
the fraction of the volume of the bed occupied by liquid, ¢ depends on the shape of the particles 
and the turbulence of the fluid. It varies with the free space of the bed at rest ¢;. 

The law is used to calculate the distance between particles as a fraction of their diameter 
and the fluid velocity. It is then possible to explain why there is no classification in a con- 
centrated bed. 

Experiments with different coloured particles show that classification is due to varying density 
and not due to variation in the free falling velocity. 


Résumé De nouvelles expériences et une étude bibliographique approfondie confirment qu'il 
y a intérét A mettre la relation entre la vitesse et la concentration du lit sous la forme u = V ¢’. 
Lexposant o dépend a la fois de la forme des grains et de la turbulence de écoulement. 
Il varie parallélement a ¢;, fraction de vide dans le lit au repos. 

En se basant sur la loi w/T — «° on ealeule les variations des distances entre grains en 
fonction du diamétre des grains et de la vitesse du fluide. On peut alors expliquer pourquoi 
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il n'y a pas de classement par difference de granulométrie dans un lit concentré. 
Des essais de fluidisation avec des solides de couleurs différentes montrent d’autre part 
que le classement se fait d’aprés les densités moyennes des suspensions et non d’aprés les vitesses 


de chute libre des grains. 


Le développement de la catalyse fluide a attiré 
attention des chercheurs sur le mécanisme de 
la mise en suspension d'une poudre par un 
courant de fluide ascendant, et sur les propriétés 
de ces suspensions ; la présente note se limitera 
a Tétude de quelques aspects importants de la 
fluidisation dans les liquides. 

Le premier probléme que nous voulons aborder 
est celui de la loi d’expansion du lit en fonction 
de la vitesse du liquide. Ce probléme fondamental 
est, de loin, celui quia été le plus étudié. Cependant 
aucune solution définitive n'est encore acquise. 

Pour définir expansion, nous employons la 
variable « qui représente la fraction du volume 
total du lit qui est occupé par le liquide. 

Nous avons proposé précédemment [3] la 
relation. 

u= Ve’ 


ot u est la vitesse ascendante moyenne du liquide, 
V la vitesse de chute isolée d'un grain moyen et 


o un nombre caractéristique des propriétés de 
la poudre, 

Nous pouvons ajouter a présent que cette 
relation est bien vérifiée par tous les résultats 
expérimentaux qui ont été publiés a notre 
connaissance [1 & 11]. Elle a le grand avantage 
détre simple et d’avoir une représentation 
linéaire dans un diagramme logarithmique. Elle 
permet la comparaison facile du comportement 
des différentes poudres et la mise en évidence 
d’anomalies dans le comportement de certaines 
de celles-ci. Elle permet également de préciser 
effet des différentes caractéristiques de la 
poudre, 

En effet, en plus des variables qui inter- 
viennent dans le calcul de V, c’est-a-dire le 
diamétre des grains et leur densité ainsi que la 
densité et la viscosité du liquide, il faut encore 
tenir compte de plusieurs facteurs qui ont une 
influence sur : 

Ces facteurs qui sont en ordre principal la 
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forme et la dimension des grains sont difliciles 
& mesurer avec précision, et nous devons nous 
contenter dindications qualitatives. 

Nous dirons que l’exposant o varie, en général, 
de 2 & 7 pour atteindre dans certains cas la 
valeur de 10. 

Les petites valeurs de o correspondent a des 
grains gros et de forme arrondie ou sphérique, 
les grandes valeurs & une granulométrie fine et 
a des formes anguleuses et irréguliéres. 

Il est intéressant de rapprocher les valeurs 
de o des propriétés de la poudre en tas que l'on 
peut caractériser en premiére approximation par 
la valeur ¢, de la fraction du volume laissé entre 
les grains au repos. 

On constate que les grains gros et sphériques 
se tassent plus fort que les grains fins et anguleux ; 
a de faibles valeurs de o correspondent done de 
faibles valeurs de toutefois une relation 
quantitative est diflicile a établir, car des facteurs 
tels que la forme des grains ou le tassement 
d'une poudre sont délicats & mesurer. 

Nous nous sommes aussi attachés & une autre 
question dans laquelle on commet souvent, 
pensons nous, de graves erreurs. C'est le probléme 
du classement qui se produit dans un lit fluidisé 
ou dans une suspension en train de décanter. 

On sait depuis longtemps que pour obtenir 
un bon classement de substances pulvérulentes 
par différence des vitesses de sédimentation il 
est nécessaire de travailler a faible concentration 
(1%, par exemple). 

On a remarqué que le classement ne se faisait 
pas bien avec des suspensions plus concentrées. 

D'autre part, dés les premiéres études de fluid- 
isation, on a été frappé par la netteté du niveau 
supérieur de lits a granulométrie trés étendue ; on 
a remarqué que cette netteté était plus grande aux 
fortes concentrations ct on a attribué le phéno- 
méne aux nombreux choes entre grains de différen- 
tes grosseurs qui auraient pour résultat d’égaliser 
leurs vitesses. 

On a toujours admis que le classement se 
faisait ou tendait a se faire selon les vitesses de 
chute libre, c’est-a-dire que les grains a grande 
vitesse de chute se placaient dans le bas, et 
ceux & faible vitesse dans le haut. 

Nous pensons qu’en se basant sur la forme 


simple de la loi qui unit vitesse et concentration : 
u = Ve’ (1) on peut donner un schéma des 
phénomeénes qui indique pourquoi la classification 
ne se fait pas dans un lit concentré. 

Pour arriver 4 ce résultat nous faisons quelques 
hypothéses simplificatrices qui nous permettront 
de calculer la distance entre les grains en fonction 
de la vitesse du liquide et de la grosseur des 
grains. 

(1) Nous supposons que nous sommes dans 
le domaine de Strokes, c’est-A-dire que la vitesse 
de chute des grains isolés V est proportionnelle 
au carré du diamétre des grains : V = k d? (2). 

(2) On considére grains sphériques 
répartis réguliérement dans l’espace, aux noeuds 
d'un réseau régulier. 

(3) On considére encore que expansion se 
fait par augmentation simultanée des distances 
entre grains sans que les positions relatives 
soient modifices. 

On en déduit que la fraction du volume 


3 
ccoupé par le solide (1 —«<) = K (3) ot D 


est la distance entre les centres de grains voisins, 
et d le diamétre des grains ; K dépend du mode 
d’empilement et vaut (1 — c’est-a-dire la 
fraction du volume occupé par le solide au 
repos, quand les grains voisins se touchent 
(d = D). 

En combinant les relations (1), (2), (3) on 


trouve : 
u=—kd (5) | 


relation entre u, d et D que nous avons représentée 
a la Fig. 1 dans le cas particulier o = 4 €, = 0-4, 
et pour quelques valeurs particuli¢res de wu. 

Ces courbes présentent un minimum qui 
permet de séparer deux zones; dans la partie 
de droite oi le rapport d/D est grand le lit 
est concentré et les courbes sont ascendantes, 
alors que dans la partie de gauche ott d/D est 
faible elles sont descendantes. 

Considérons d’abord la partie de droite, pour 
une vitesse u, = 0-1 des particules de 75y sont 
en équilibre & une distance de 95x. 

Imaginons parmi les autres, une particule 
plus petite B, pour qu’ elle reste a la méme 
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20 30 40 30 60 70 80 90 100 110 

Fig. 1. Relation entre le diamétre de grains d, la distance 

entre les grains D, et la vitesse moyenne du liquide u. 
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distance de 95 de ses voisines il faudrait une 
vitesse u, plus grande que u ,. 

Cela veut dire que si la vitesse du liquide est 
u, cette particule va descendre avec une vitesse 
Up — U4, ce faisant elle se raproche de ses voisines 
inférieures, D diminue et on atteint un nouvel 
équilibre représenté par le point C. 

On voit done qu'une particule plus petite au 
lieu d’étre entrainée par le courant et de monter 
comme cela se passe aux faibles concentrations 
tend dans ce domaine aA se rapprocher de la 
couche inférieure; inversement une particule 
plus grosse G trouve sa position d’équilibre en 
H & une distance plus grande aprés avoir monté 
puisque u, est plus grand que uy. 

On peut illustrer cette explication mathé- 
matique et qui n'est évidemment que schématique 
par l'image suivante. 

Considérons quelques grains voisins, supposons 
que le grain central devienne brusquement plus 
petit, il y a dans son voisinage une augmentation 
de la section libre qui provoque une diminution 
de la vitesse locale réelle du fluide et done une 
diminution de la force d’entrainement. II se 
peut que cette diminution soit plus que suffisante 
pour compenser la diminution du poids apparent 
de la particule; dans ce cas la particule descendra 
et viendra trouver une nouvelle position d équilibre 
plus proche de ses voisines inférieures. 


Dans la partie gauche de la figure, aux faibles 


concentrations, les courbes sont descendantes, 
les phénoménes sont inverses, les grains les plus 
gros ont tendance A descendre, les grains les plus 
fins & monter; Ceux-ci peuvent méme étre 
entrainés en dehors du lit, qui se dissociera alors 
en plusieurs zones. 

Les points minimum des courbes correspondent 
aux valeurs limites de Texpansion au-dela 
desquelles il n'est plus possible d’avoir d'une 
facon stable un lit homogéne avec une limite 
supérieure bien nette. 

On peut montrer qu'on a 


Cmax = (1 «) pour =te = O04 


on trouve ¢,,., = 42 ce qui est en bon accord 
avec l’expérience. 

Le mécanisme exposé ci-dessus n'est valable 
que pour décrire le classement de grains tous de 
méme densité. 

L’étude expérimentale que nous avons entre- 
prise pour en vérifier les résultats nous a cepen- 
dant permis de tirer quelques conclusions con- 
cernant le cas de deux matiéres différentes, 
ayant des densités différentes. 

Nous avons fluidisé dans l'eau des mélanges de 
poudres de couleurs différentes, mais de dimensions 
bien limitées et dont nous connaissions les 
vitesses de chute libre; de cette fagon nous 
pouvons suivre aisément le classement. 

Nous pouvons décrire les phénoménes observés 
de la facon suivante : 

ler cas. On considére un mélange de deux 
poudres de méme densité mais de granulométries 
différentes (nous avons employé deux fractions 
de sable de quartz, classées par lévigation et dont 
une a été colorée en bleu). 

Il_n’y a classement que si la différence des 
diamétres moyens est suflisante. (telle que 
V, = 2 -V,). Pour des différences de granulo- 
métrie plus faibles, et aux fortes concentrations 
il n’y a aucun classement, le lit reste homogéne 
aussi longtemps que dure l’expérience. 

Nous avons la une confirmation du processus 
expliqué ci-dessus : il n’y a pas de classement. 

Dans le 2¢ cas, on considére un mélange de 
deux poudres de densités différentes (par exemple 
du sable et du charbon, ou bien du sable et du 
ferro silicium). 
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Il y a trois possibilités : 

(a) La vitesse de chute libre du solide le plus 
lourd est la plus grande. On observe alors que le 
classement se fait rapidement, le solide le plus 
dense et A grande vitesse de chute occupe le 
bas du lit, le solide le plus léger et a plus faible 
vitesse de chute occupe la partie supérieure, 
aprés un délai de quelques minutes la limite 
entre les deux zones devient définitivement bien 
nette. 

Ce classement a licu méme si les vitesses de 
chute libre sont trés proches lune de l'autre. 

(b) Une deuxiéme possibilité est que les 
vitesses de chute libre soient identiques. les 
diagrammes « = f (u) des deux poudres sont 
confondus. 

Dans ce cas on observe qu'il y a un classement : 
le solide le plus dense va oceuper la partie inférieure 
du lit. On ne peut pas expliquer ce classement 
par une différence entre les vitesses de chute des 
grains. 

(c) Dans le troisiéme cas, ot la vitesse de 
chute libre des grains denses est plus petite que 
celle des grains légers, le classement se produit 
encore: les grains denses a faible vitesse de 
chute occupent la partie inférieure du lit. Tout 
se passe comme si les deux poudres formaient 
des phases séparées et se comportaient comme des 
liquides non miscibles. 

On peut expliquer ce phénoméne par le méme 
processus que précédemment. 

Les grains de faible densité ont une plus grande 
vitesse de chute libre, ils sont done beaucoup 
plus gros que les autres, il y a dans leur voisinage 
une forte réduction de la section libre offerte au 
liquide qui provoque une augmentation locale 
de la vitesse, suffisante pour les entrainer vers le 
haut, malgré leur plus grande vitesse de chute 
libre. 

D’autre part, un classement en sens inverse, 
qui conduirait la poudre la plus légére dans 
le bas du lit est impossible puisqu’il donnerait 
lieu & une situation instable avec une phase plus 
légére en dessous d'une phase plus lourde. 

Une derniére expérience met bien en évidence 
limportance de cette question de stabilité. 

On mélange deux poudres dont les vitesses de 
chute sont suffisemment différentes pour que 
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Fig. 2. Représentation de la fluidisation simultanée de 


deux solides. 


les domaines de fluidisation ne se superposent 
pas : (Fig. 2). 

La poudre la plus dense (1) a une plus petite 
vitesse de chute. 

Pour une vitesse u,, la poudre 2 n’est pas encore 
fluidisée, elle se sépare néanmoins de la poudre 
1 qui, elle, est bien fluidisée & cette vitesse, mais 
chose assez remarquable elle occupe la partie 
supérieure du lit. 

Les grains sont trop gros pour étre soulevés 
individuellement et mis en mouvement par le 
liquide & cette vitesse, mais l’ensemble du lit 
encore compact vient flotter au-dessus du lit 
fluidisé formé par les grains plus denses. 

On voit done que le facteur primordial dans 
cette étude du classement des lits fluidisés est la 
densité moyenne des phases fluidisées formées par 
chacun des constituants du mélange pris séparé- 
ment. 

Si nous augmentons la vitesse du liquide, le 
lit de solide 1 continue a gonfler, a s’étendre, 
sa densité moyenne diminue, pour une vitesse 
u, cette densité devient plus faible que celle du 
lit 2 qui ne varie pas, on observe alors un renverse- 
ment, aprés une période ott les deux poudres sont 
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mélangées un nouveau classement se fait dans 
lequel la poudre 2 occupe la partie inférieure 
et la poudre 1 la partie supérieure du lit. 


ne serait-ce que comme un outil de travail 
constituant une premiére approximation. 
Ce travail a été effectué au Service de Chimie 


Industrielle de [Université Libre de Bruxelles, 
sous la Direction de Monsieur le Professeur 
De Keyser, avee laide d'une subvention du 
Fonds Tassel. 


Le calcul montre que ce renversement se 
produit bien & la vitesse ot: la densité des deux 
lits devient égale. 

Un mélange de deux poudres peut done selon 
la vitesse du liquide se classer de deux facons 
diamétralement opposées, cela prouve que le 
classement dépend des propriétés globales des d = diamétre d'un grain. 
suspensions et non seulement des propriétés |  D distance entre deux grains. 
particulieres des grains. ¢ = expansion du lit, c’est-d-dire rapport de la hauteur 

En résumé nous dirons que le classement meine 

Waprés les granulométries doit se faire dans 1, = hauteur du lit au repos. 
une phase diluée ; dans un lit concentré le classe- u = vitesse du liquide au-dessus du lit fluidisé. 
ment se fait selon les densités. I’ = Vitesse de chute libre des grains isolés dans l'eau. 
liquide. 
en dehors de son intérét propre permet dillustrer 
en quelque sorte lintérét qu’ il y a & adopter la 


NOTATION 


fraction du volume occeupé par le liquide dans 


le lit au repos. 1954 
forme simple de la loi d’expansion u = Ve’ o ~ exposant de ¢ dans la relation empirique a = Vee 
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Summary—In the evaluation of thermodynamic measurements of mixtures a development 
of the excess Gibbs free energy, AG, in orthogonal functions has proved advantageous. Since 


measurements do not give a smooth curve but a set of discontinuous values, a method appropriate 
for the evaluation of equidistant discontinuous measurements has been developed, using Verrin’s 
discontinuous orthogonal polynomials, The necessary formulae and tables are listed in this paper. 


Résumé— Pour lestimation des propriétés thermodynamiques des solutions, il est intéressant 
de développer en séries orthogonales la fonction AGE, énergie libre d’excés de Grass (ou enthalpie 


libre d’excés). 


Comme les expériences ne donnent pas une courbe réguliére mais une série de valeurs discon- 
tinues, il a fallu développer une méthode convenable pour lévaluation de mesures discontinues 


équidistantes, en recourant aux fonctions polynomiales discontinues de Verrin. 


et tables nécessaires dans le mémoire. 


The method generally adopted in the evaluation 
of equilibrium measurements of binary liquid 
mixtures [1] consists in a power series develop- 
ment of the molar excess (Gibbs) free energy 
This series is related to the so-called 


AGE, 
Margules [2] series. Its most appropriate form, 
probably, is the symmetrical development, 


introduced by Repuicu and Kisrer [3]. 


AGE — B &, &, a=2f, —1=§, —,. (1) 
The determination of the constants a, becomes 
particularly simple if special values of &€ are 
preferred in the evaluation. This advantage, at 
the same time, implies a weakness of the method : 
the constants are not determined in such a way 
as to give, on the average, best agreement with 
all experimental points. 
Our aim has been to find a method of repre- 
senting the results which 
(a) gives on the average optimum agreement 
with all experimental points, and 
(b) furnishes a final determination of the con- 
stants entering the development. 
By (b) the following is meant: if the coeffi- 
cients of the first n members of the development 
have been determined, and if it is found out that 
this development is unsatisfactory then we try a 


Formules et 


the first » constants, however, need not be deter- 
mined anew, but retain their original values. A 
development in orthogonal functions complies 
with our requirements, as, for instance, a Fourier 
series development, which, however, is not the 
suitable solution for our problem. Considering 
that the first two terms in the development (1) 
correspond to P, and P,, the first two Legendre 
polynomials (spherical harmonics) [4] of 
w= 2£—1, varying from —1 to +1, the 
obvious development is 
AGE 

BE & 
Here the Legendre polynomials of the first kind 
are defined by 


(2) = — UY, (3) 
the first 6 polynomials being 

P,(v) =1 

P, («) =2 

P, (x) = 3 (82? — 1) (+) 

P, (x) = 4 (52° — 82) 

P, (x) = § (85 a+ — 302? + 3) 

P, (x) = § (682° — 702° + 152) 
The coeflicients A, are given by the relation 

+1 
de 


development with m > n members of the series ; 
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which corresponds to the analogous equation for 
a Fourier development and is easily understood. 

This development is satisfactory in every 
respect, if a continuous curve (for instance by 
graphical interpolation) for the experimental 
values is available. Especially, it allows of an 
unambiguous, straightforward determination of 
the coeflicients. 

If, however, the experimental points are not 
sufliciently closely spaced, systematic errors still 
may be introduced by the arbitrariness left in 
drawing a smooth interpolating curve. Conse- 
quently, we were looking for a procedure offering 
the advantages of a development in orthogonal 
functions, but making use only of a discrete set 
of measured values. At this point the advice of 
Professor A. Warner, Darmstadt, was of 
decisive help to us. He drew our attention to 
Verrin’s[5] discontinuous orthogonal polynomials, 
by which the problem is completely solved. 

The authors are indebted to Professor van 
KREVELEN for calling their attention to the fact 
that analogous orthogonal polynomials are used 
in statistical theory (R. A. Fisner, “ Statistical 
methods for Research Workers "; R. A. Fisner 
and F. Yates “ Statistical Tables,” Oliver & 
Boyd, London). The tables as given in this 
paper however are more convenient for the 
present purpose. 

These polynomials V,(# w) are a generaliza- 
tion of Legendre’s polynomials, valid for a dis- 
crete set of equidistant measurements, defined in 
such a way that 


lim = P, (x) 


w—>0 


(6) 


where P (x) is the Legendre polynomial, w being 
the length between equidistant experimental 
points in Fig. 1. We now attempt a development 
for our function f (2) (equation 2) 


(7) 


fla)= Za, V, (aw). 


The coefficients a,, corresponding to the coefli- 
cients A, of the former development in Legendre 
polynomials, are given now by sums (instead of 
the former integrals) 


#2 +e 


Discontinuous 
orthogonal 
polynomials 
w= Ar=02 


First derivatives of 
the discontinuous 
orthogonal 
polynomials, 

w= Ar = 0-2, 


Fig. 1. 


2-(A—l)e 


where N = 1,0,» =a w = N-a, 2N + 1 being 
the number of equidistant experimental points 
(including the two points corresponding to the 
two pure components). » and N are introduced 
only as variables for the summation. The first 
Vettin polynomials are 


24+ 


+ w)) 


1 
(10.2% — (6 + 6 w—2w*)r) 
1 
35a* —(30 4+ 30w 


+ (8 +6 w —38w? —6 (9) 


If the functions V,(2/w) and the factors f,, are 
tabulated (cf. the tables, App. I). 
(2v+1)0 2—(A—1)a 

Awo2+(A 
then the evaluation of the coeflicients a, requires 
a summation of products only : 2 f(y) V, (»). 
which is quickly done by means of any computing 
machine. The necessary tables are very short, 
because values for a discrete set of arguments only 


are needed, 


2+ 


Fig. 2. 


The only drawback of this method, still remain- 
ing, consists in the fact that equidistant experi- 
mental points are required. This, however, is not 
of a serious nature. If one carries out measure- 
ments, it is easy to have two or more points in 
the vicinity of, say, 10, 20 . . . . 80, 90 mol %, 
and to interpolate without appreciable error for 
the exact values 10, 20, . . . . 80,90 mol °. In 
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the examples evaluated by us, using published 
results, we had to employ a graphical interpolation 
(of AG®/B &,), ef. below, p. 22. The remaining 
small deviations between experimental and caleu- 
lated values (of usually less than 1°%,) are probably 
due to this interpolation. 

The constant B in equation (1) or (2) is arbi- 
trary. It is convenient to use B 4 AGF, 
because then the function f (7) assumes values in 
the vicinity of unity, ef. the typical curves in 
Fig. 2. Of these curves 1 and 3a are easily 
approximated by means of not more than 5 
Vettin polynomials (up to V,), the approximation 
for 2 is somewhat less good, more difficult is that 
of 3b and e. 

The Gibbs-Duhem equation imposes no res- 
trictions upon the function f(#). Therefore it is 
legitimate to replace equation (2) by 


BE, 
where g(x) is a non-vanishing arbitrary simple 
function of x, chosen in such a way as to give 
f(x) a shape which is conveniently approximated 
by not too many Legendre or Vettin polynomials. 
This procedure is quicker than the use of more 
polynomials for an improvement of the approxi- 
mation. 

The same developments may be used for the 
evaluation of measurements of total pressures of 
binary systems where measurements of partial 
pressures are lacking. A decomposition of the 
total-pressure curve into two curves for the 
partial pressures is usually done by an integration 
of the Duhem-Margules equation. It is, however, 
much simpler and more straightforward to use a 
development for AG®, consistent with the Gibbs- 
Duhem equation, and to determine the constants 
in such a way as to give best agreement on the 
average with measured total pressures.* This 
method is applicable with Vettin polynomials, 
with Legendre polynomials or with a simple 
power series development for f(x). Details will 
be published elsewhere. 


(2) (11) 


~ * A special case of this procedure has been applied by 
A. Musi [6], who used the first two or three members of 
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the Margules series. 
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APPENDIX I Table 1. V, for w = 0-2; v = 1, 2, 3, 4, 5. 
. , | j 
fa) | The V, (w/w) for w = 0-2 and | slats 
V, = d/dx(V,) 
% = o | | 1] e 
o2 | 15 |-3/s |-7/15| 4/8 
10 0-4 2/5 | | 23/30, | 4/3 
) 06 3/5 | -2/% -11/15) -1 | -1/8 
9 5 3 4/5  +4+2/5 | -1/5 —1} 
Vs = — 7122) (42 0-89) | | | | 
V5 = — 352" Table 2. V', for w = 0-2; = 2, 3, 4, 5. 
+ 4032) 
= (12625 — 15962" = = (6300*— 478-827 | 2 | 3 4 | 
88-4384 2) 88-4384) | | 
0 0 — 39/86 0 3003 /378 
(b) the same for w 0-1 o2 2/3 37/18 115/36 1057 /252 
4/3 20/36 — 85/18 — 476/63 
06 2 23/18 — 35/12 |—2723/252 
- O8 8/3 151/36 35/9 — 259/126 
0 0 _—2-47223 | 0 
658 
"sin *) o2 066666 —205555 
300 O+ —0-80555 —472222 — 755555 
Vg = — + 06 200000 127777 —2-91665 18-8055 
O8 2-66665 419444 — 205555 
= (12645 — 151-929 + 35-2024.) 1 3-33333 17 -B611T 39-10464 
S814 
(462.2° — 678-32" + 242-6002? — 12-3552) 
6" 8721 Table 3. The factors f,,, for w = 0-2. 
7= (171627 — 2956-8 2° + 1428-432 — 168-093 2) » | Say 
In the following the polynomials are tabulated ° 1/11 = 0-0900080 
for 0 < + 1, the values for — 1 < < Oare 
obtained by considering that V,, is an even, 3 30/143 = 0-20979 
V,,,4, an odd function of x. 4 18/143 = 0-12588 


Table 4. for w = 0-1; = 1, 2, 3, 4, 5, 6, 7. 
2 3 4 | 5 6 


—O0-57895 0 0-61300 0 | — 0-88545 0 
— 0-56316 18947 O-55728 036223 — 0-71641 — 0-78885 


7 


O-2 —O-51579 — 036140 O-39732 0-63055 — 026563 — 118328 
0-3 — 0-43684 — 049825 015480 0-72730 0-30975 — 095572 
O-4 — 032632 — 058246 13414 0-60733 O-78762 — 0-16810 
OS —0-18421 — 0-59649 — 0-41899 027425 094041 O-79414 


— 0-01053 — 052281 — 0630038 — 020330 061858 132077 
O7 0-19474 — 034386 —O-70O175 — O-67544 — 014210 0-85267 
O-8 043158 —0-04211 — 052632 —O-89474 — 0°98946 — 0-62100 
0-70000 0.40000 0 —0-50000 —1:10000, — 1-79992 
1 | 1 1 1 1 1 1 
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Table 5. Factors f,, for w = 0-1; v = 0,1, 2, 3, 4, 5, 6, 7. 


Say 


0-047619 
0-12987 
0-17881 
0-18775 
0-16414 
0-12346 
0-081059 
0-046765 


S 


Sometimes it is desirable to rewrite a series of 
Vettin polynomials as a simple power series. We 
give the formulae (a) for n — 4 and w = 0-2 


i= ZA,2’ 


A,= 


(b) forn = 7, and w = 0-1 


= dy — 057895 dy + 0-GO944 ay — 088545 ag 
= — 192399 a, + 378426 a, — 860468 a, 
157896 dg — 560025 ay + 17-B8682 ay 

= 2924 a, — 16-32925 a, 73-12133 a, 

= 5985 ag — 48-61105 ay 

13-545 a, — 151-35859 a, 

= 33-10969 ag 

A, = 8784204 a, 


It is seen that a change in the even coeflicients 
ay, affects the even coetlicients Ay, only, and vice 
versa, 


Aprenpbix Il 


Detailed evaluation of an example. The system 
ethyl aleohol-toluene was chosen and measure- 
ments of Wiese and Krerscumer [7] were 
evaluated. The results obtained by the authors 
are tabulated below, cf. Fig. 3. 


-6o 
° 


© MEASURED VALUES 
@ CALCULATED VALUES 


Fig. 3. Thermodynamic functions of the system ethyl 


aleohol — toluene at 35°C, experimental values after 
Krerscumer and Wiese. 


Table 1. Experimental values for ethyl alcohol-toluene, after 


Krerscumer and Wiese. 


AGE f(z) | TASE 
557 | 55° 35° 


AGE f(x) 
35° | 33° 
00439 «63 
O-115T 143 
236 
04043 280 
281 
06282 256 
O-7186 219 
O-8423 «143 
09163 83 
0-9635 37 


1-204 
1-204 
1-088 

1-003 
0-999 
0-044 
0-933 
0-928 

| 0-933 

0-906 

| 


—25 


—12 


Act was calculated for 
BE, 


35°C with B = 1160, for 55°C with B= 1164. 
From the curve for f(«) values for € = 0, 0-1, 
0-2... . 0-8, 0-9, 1-0 were interpolated, Table 2. 

Now the coeflicients a, of equation (8) using 
the notation of equation (10) must be deter- 
mined. For this purpose the sums must be 


The function f (7) = 


a 
| 200 
| 
| 
| | 
£ 
100 
le 
| 
« 
4 Ag = — + % 
89 
=a, ——a 
5 625 
Ay = 
2 32 76 
125 
36 3 
625 
Bed 
150 68 211 
232 «(1-064 36 | 272 
281 1002 — 278 
282 0-999 — 4 277 
258 0-948 —35 221 
221 0-938 —441 178 
147 0949 —35 108 
84 OOF] 58 
38 (0-927) | 25 
21 
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computed, for instance for the determination of 
a, (with w = 0-2) 
a = V, (2) 
the sum according to (8) 
2 f (2) - V, (x) =f(—1) (—1) 
+J(— 0-8) ( — 0-8) 
+ f(08) (08) 
+fQ) V0) (11,2) 
Table 2. Interpolated values for f (x), from values of Table 1. 
| | | | 
| 


| 1381 


(11,1) 


1-218 
1-044 
1-002 
0-971 
0-950 
0-936 
0-929 
0-919 
0-909 


—3 
1 
2 
3 
4 
5 


— 


| 


— 


Fig. 4. Example of approximation by means of discon- 
tinuous orthogonal polynomials ; AG for ethyl alcohol- 


These sums and the corresponding coefficients a, 
for 35°C and 55°C are listed in Table 3, for 
v = 0,1, 2,3, 4. 


Activities of i 
ethyl alcohol- 
toluene, 35°C, 
calculated by 
means of AGE 

= BE, & fa,V,. 


Activity coeffici- 
ents of ethyl 
alcohol-toluene, 
35°C calculated by 
means of AGE 

= BE 


i7-- 
Fig. 5. 


In the following table the results obtained for 
the first, second, third and fourth approximations 
(up to vy = 1, 2, 3, 4) are compared with experi- 
mental values f(x), and for the fourth approxima- 


ce’, is given in the last 


tion the ratio «, = Tia) 


column. (¢€, — 1) 100 is the percentage deviation, 
and it is seen that the maximum at 35°C is 
0.35%, at 55°C 0.83%. 


| 
0 1321 | -1 
1 1168 | 
o2 | 1087 06 
Os 1038-04 | VOL. 
| «1002 —02 3 
0-975 0-0 | 
os 0-957 o2 logy, log yz 1954 
0-045 O06 
oo 
| 09300 10 
° 
of + | | 
toluene at 35°C. ee 
Table 3 
| a, Lf a). V, “, | 
| 35° 33° | 55° | 55° AT 
o 11-376 | 10342 | 10284 
| —0-8776 | —0-1995 | ~0-70382 —0-1598 1-982 10-8 
2 | 0-3979 0-1043 O-3385 0-08878 —7-78.10°* 
s | ~0-1673 —0-0351 —0-1520 — 0-03189 | 161.104 
| 0-04933 0-00621 0-06081 0-00766 0-73. 10°5 
22 
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The «,’s, equal to the ratio of calculated sums 
and f(x) as obtained by interpolation between 
values derived from experiment, show the 
rapidity of convergence of the method. In order 
to show this more figures have been retained than 
are justified by the accuracy of the values AG® 
listed in Table 1. It should be kept in mind, 
however, that any function f (7) or AG® may be 
approximated in this way. For a comparison 
with experiment and for a test of the thermo- 
dynamic consistency of the experiments it is 
necessary to calculate backwards the individual 
P, (or Ap*) or the total pressure P. In this 


43-00 
1-02 42.00 
rot oo 
0-99} 
onal > 
°-98 anes ° 17-200 & 
13.00 * 


Fig. 6. Deviations from experimental values for several 
approximations by means of discontinuous orthogonal 
polynomials AGE for ethyl aleohol-toluene at 35°C. 


Table 4. 
| | 
0 0 

= 1-381 | 1-2336 1-3378 1-3731 1-3838 0-9988 
1-218 1-1936 1-2355 12285 12142 1-0035 
~06 1-1537 11469 0-9982 
1-044 | 1-1140 1-0722 10458 1.0407 10003 
1-022 10741 10115 0-9951 1-0046 0-9973 
0-0 0-971 1-0342 09646 09646 060-9700 0-9999 35°C 
0-2 0950 0-9317 09481 1-0023 
0-936 09544  0-9127 0-9396 0-9332 10027 
0-6 0-929 0-9145 0-9076 0 0-9284 09979 
0-8 0-919 0-87.46 0-9164 0-9234 0-9252 
10 0909 @ 0-9391 0-9040 0-9057 10013 

| | | | 

| 2 3 | 4 
z= Zz Zz % 

1168 1-1917 1/1854 1-1777 1-0083 
| 1087 11183 10950 00083 
1-038 0928 10568 BT | 1-0004 
1002 10603 10071 0-9922 09973 ©9953 
oo 0-975 1-0284 09692 09692 09768 10019 55°C 
| 0-957 0-9964 09431 9580 09631 
0-949 09644 06-9289 0-9534 09521 11-0033 
06 09325 0-266 0 0 0-9971 
os 09005 0942 0-948 0-99-44 
10 | 0-930 | 0-8686 0-9573 | 0-9254 | 0-9331 | 1-0033 


procedure in addition to the absolute values of 
J (2), its derivatives, too, appear. The accuracy 


obtained may be judged from Table 4a. 


It should be considered however that deviations 
may partly be due to inaccuracies involved in the 
interpolation. 
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Table 4a. 


Preate 
77-63 
84-90 

102-01 

107-93 

113-49 

114-99 

117-63 

118-22 

16-11 

107-34 


0-0330 
0-0468 
O-1214 
02079 
0-3620 
0-4160 
59300 
0-7263 
O-8519 
0-9701 


117-00 
118-57 
116-56 
107-64 


Pons Preate 


+ 144 
+ 0-08 
+100 
+ O35 
+ O27 
+ O35 
+ O45 
+ 0-30 


"cate £1 “obs ‘cate 


0-0078 
0-0062 
— 0-0048 
— 0-0016 
00027 
—O-0110 
—0-0079 
00-0015 
—0-0012 
0-0017 


O-4138 
0-4687 
0O-5710 
0-6030 
0-6373 
0-6809 
O-7149 
O-7836 


O4749 
0-5662 
0-6014 
00-6546 
0-6384 
0-6730 
O-7164 
O-7848 
O-9318 


£,; = mole fraction ethanol in liquid 
£, = mole fraction ethanol in vapour 
P = total pressure, Torr, at 35°C 


cale = 


In addition the functions 


0, = (11,3) 


T. ASE — £,& TBS 0, V.; 
0 

and 

(11,4) 


a, + To, 


AH = 
0 


d (4, A, 
(3) --* (11,5) 
were computed. Here B was considered as con- 
stant, which is legitimate because it changes less 
than 0-4%, from 35 to 55°C. The last equation 
(11,5) allows of a calculation of the a,’s for other 
temperatures, according to 


The equations thus obtained are 


(11,6) 


AGE = 1160- —0-1995 40-1048, 
a” 0-0351 + 0-00621 


AGE = — + 0-08878 
55° — 0-03189 + 0-00766 (11,7) 
— ASE = 0-1160 419-821, —7-78 Vy 
+ 161 Vy + 0-073 Vy} 
— 00145 + 0-0039 V,} 


obs = observed values 
calculated values with 


4 
= Za,V, 
0 


The numerical values for 35°C are listed in 


Table 5. 


Table 5. Thermodynamic functions for the system cthyl 


alcohol-toluenc, at 35°C. 


The complete evaluation of this system, by 
means of a normal calculating machine, took 
about 4-5 hours. Included in this time is the 
repetition of all important steps of the calculation, 
for control. 


ACKNOWLEDGMENT 
The authors are indebted to Professor A. 
Wattruer, Darmstadt, Institut fiir Praktische 
Mathematik, for calling their attention to the 
discontinuous orthogonal functions. 


24 


Pos 
79-38 0-4216 
86-34 
102-09 
108-93 
126 
115-34 
VOL. 
3 
1954 
AGE TASE Mu 
0 0 0 0 
o-1 127-6 69-6 197-2 
0-2 206-9 75-3 282-2 
0-3 2344 48-8 303-2 
278-2 120 290-2 
0-5 281-5 ~210 260-5 
0-6 265-1 222-1 
0-7 228-6 —31-7 176-9 
0-8 172-1 —~472 124-9 
oo 95-8 — 65-4 
0 0 0 0 
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NoTaTIOoN 


AGE = molar excess Gibbs free energy, cal /mol 
B = a constant of dimension, cal/mol 
£,. &g = mole fractions of components 1 or 2 respectively 
# = symmetrical variable, # = 2£, — 1 
a, 4, = constants 
P,, = Legendre polynomial of kth order 
V, = Vettin’s polynomial of vth order 
ASE = molar excess entropy, cal/mol °K 
AH = molar heat of mixing, cal/mol 
o,, A, = constants 
Note added in proof : 
A complete survey of the problem of » least squares filting 
of data by means of polynomials” has been given by 
Birce, R. T. and Weinperc, I. W.; Rev. Mod. Phys. 
1947 19 298. 
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The distribution of residence-times in an industrial fluidised reactor 


P. V. Danckwerrts, J. W. Jenkins and G. Place 


Department of Chemical Engineering, Tennis Court Road, Cambridge, England 
( Received 8 February, 1954) 


Summary—The distribution of residence-times of the gas flowing through a fluidised catalyst 
regenerator has been determined by using helium as a tracer. Only about 5% of the gas spends 
less than half, and about 5% more than twice, the mean residence-time in the regenerator. It 
appears that the flow through the fluidised bed itself is much closer to piston-flow than to complete 
mixing, the hold-back being less than 0-1. 


Résumé—En utilisant lhélium comme indicateur, les auteurs ont étudié la distribution des 
“ durées de séjour d'un gaz qui s‘écoule au travers d'un régénérateur de catalyseur fluidisé. 
Tandis que 5% de la masse totale du gaz présente une “ durée de séjour ™ inférieure a la moitié 
de la valeur moyenne, il y en a aussi 5% dont la “ durée de séjour * dépasse le double de cette 
valeur. I] semble qu’au travers d'un lit fluidisé, 'écoulement gazeux se rapproche beaucoup 
plus d’une veine cylindrique rigide que d'un régime turbulent, la “ rétension " reste inférieure 


a 0-1. 


When gas flows through a fluidised bed of solid 
particles, there is a considerable spread of 
residence-times ; different molecules of entering 
gas spend different lengths of time in passage, so 
that part of the gas spends less, and part more 
than the average length of time in the bed. This 
spread of residence-times, which is caused by 
back-mixing, or departure from piston-flow, is 
one of the most important characteristics of a 
fluidised reactor, as it will influence the nature 
and rate of the reaction. As yet, however, there 
is very little published information on the subject. 

GILLILAND and his co-workers [3] have meas- 
ured the distribution of residence-times (and 
related characteristics) in model fluidised beds. 
However, there is no certainty that these repro- 
duce the behaviour of the large-scale beds used 
industrially, For instance, in small beds the solid 
material can be seen to follow a circulatory path, 
with dimensions determined by those of the con- 
fining vessel, and this circulation of the solid 
helps to bring about back-mixing of the gas [3]. 
However, it is not possible to predict the way in 
which this motion will develop as the absolute 
dimensions of the bed are increased. It would be 
exceedingly difficult, perhaps impossible, to 
ensure rigorous dynamical similarity between 
fluidised beds whose linear dimensions are in 
the proportion of, say, 120 to 1 (40 ft to 4 in.). 


Such similarity was not attempted in GILLILAND’s 
work, 

Some of the difficulties involved in the use of 
models can be avoided by carrying out experi- 
ments on industrial equipment. Other difficulties 
then arise, however; for instance, the points at 
which measurements are made are often dictated 
by their accessibility, and it is not usually possible 
to cover a great range of variables. Nevertheless, 
information obtained in this way may provide a 
wholesome test of the value of data obtained in 
the laboratory, as well as being specifically 
applicable to the equipment on which the 
experiments are made, 

The only published work which provides 
information on the distribution of residence-times 
in a large fluidised bed is that of Asxrns et al [1]. 
Samples of gas were taken from a number of 
points in a cracking-catalyst regenerator. The 
uniformity of composition of these samples 
suggested that the gas in the bed was completely 
mixed ; however, there was reason to suppose 
that the samples were not truly representative of 
the gas passing the sampling-points, since pro- 
portionately less gas was drawn from the “ bub- 
bles which occur in the bed than from the 
interstitial gas flow. This phenomenon was also 
encountered by G1ILLILAND and Mason [8]. 

The experiments described here were carried 
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out on a catalyst regenerator virtually identical 
with that used by Askrns et al [2] although the 
air-rate and quantity of catalyst may have 
differed in the two cases. 

A study of the residence-times of gas in a 
blast-furnace, using radon as a_ tracer, was 
carried out by Voice [6], and has some points 
of similarity with the work described here. 


GENERAL PRINCIPLES 


The problem of defining and measuring the dis- 
tribution of residence-times in continuous-flow 
systems has been discussed elsewhere by Danck- 
werts and by other authors ([2], [3], [4], [5]). In 
outline, the method adopted here was to inject 
a dose of tracer (helium) at the entrance to the 
system, and to determine its time-distribution in 
the outflowing stream. 

Fig. 1 is a diagrammatic view of the regenerator. 
Spent catalyst from the reactor is blown up the 


stack 


10 20 30 40 ft 
SCALE 
@ 

ain HA 


Fig. 1. Fluidised catalyst regenerator. 


catalyst riser by a stream of air and enters the 
bottom cone of the regenerator. Between the 
bottom cone and the main body of the regenerator 
is a grate-like grid. The area of the slots in the 
grid which allows the passage of gas and catalyst 
is approximately 10%, of the total cross-section 
of the cylindrical part of the regenerator. Catalyst 
drawn from above the grid passes down two pipes 
to the catalyst recycle coolers (only one is indi- 
cated on the diagram); the catalyst is blown up 
through the coolers into the bottom cone by a 
stream of air. Roughly half the air enters the 
regenerator via the main riser, the remainder 
entering with the recycled cooled catalyst and 
through the pipe A. The regenerated catalyst is 
withdrawn from the space above the grid, and 
the gases pass out via the cyclones and the 
stack. 

The procedure in these experiments was to 
inject rapidly a measured amount of helium into 
the air-line to the catalyst riser, and subsequently 
to take a series of samples of the stack gas at 
measured intervals. The ratio of helium to 
nitrogen was determined in these samples. 
Helium was used as a tracer (a) because of its 
unique molecular weight, which enables it to be 
estimated easily by mass-spectrometry;  (b) 
because there is a negligible amount normally 
present in the stack gas; (c) because it would not 
be adsorbed by the catalyst. 

The amount of information which can be 
derived from such an experiment depends on the 
validity of a cumulative series of assumptions. 

(a) Assuming no gross experimental errors, and 
provided the injection is completed in a time 
much less than the mean residence-time in the 
regenerator (c. 30 secs) (a condition fulfilled in 
these experiments), the experiment gives unequi- 
vocally the distribution of residence-times of the 
helium atoms in their passage from the injection- 
point to the sampling-point. 

(b) Assuming that molecular diffusion and 
gravitational convection play a negligible part in 
determining the paths of gas-molecules in the 
regenerator, in view of the high flow-velocities 
reigning in most regions; and that the injected 
helium is uniformly and rapidly spread over the 
cross-section of the air-pipe by turbulence 
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(Re 10°); then the distribution of 
residence-times of the nitrogen which enters via 
the riser will be the same as that of the helium. 
The helium can then be regarded as * labelling ” 
the nitrogen with which it enters. 

(c) Assuming that the spread of residence-times 
between the injection-point and the bottom of the 
regenerator is very small (i.e. all the helium atoms 
enter the regenerator in a time much less than the 
mean residence-time in the regenerator) and that 
the time of passage between these points (c. 5 sec) 
can be calculated, it is possible to find the distri- 
bution of residence-times of the nitrogen in the 
regenerator itself. 

(d) Virtually all the air entering the regenerator 
does so by roughly the same path, as it enters via 
the bottom cone below the grid. The stream of 
air and catalyst entering via the riser is split by 
a distributor into a number of non-vertical 
streams. The interaction of these with the streams 
entering from the catalyst coolers and by pipe A 
no doubt gives rise to a high degree of agitation 
in the space below the grid. If it is assumed that 
this brings about perfect mixing in this space (i.e. 
entering gas molecules are rapidly and uniformly 
dispersed throughout the space) then the distri- 
bution of residence-times will be the same for 
nitrogen entering by any of the various routes. 
The experiment can then be taken to give the 
distribution of residence times of the nitrogen 
flow as a whole in its passage through the regenera- 
tor. 

The quantity of nitrogen leaving with the 
regenerated catalyst is known to be negligible in 
comparison with the total flow, so that virtually 
all nitrogen leaves through the stack. A small 
amount of fluidising air is put into the down- 
comers to the recycle coolers, and it is considered 
that this prevented any significant amount of 
helium being carried round the recycle cooling 
circuits ; the residence-times measured therefore 
refer to the regenerator exclusive of the coolers. 


EXPERIMENTAL 
Helium Injection— The helium was introduced from the 
storage cylinder into a pressure-vessel from which air had 
been displaced by filling it with water. The volume of 
gas in the vessel was about 1-3 cu. ft. and the pressure 
about 150 p.s.i.g. before and about 50 p.s.i.g. after helium 


had been released into the air-line. The pressure-vessel 
was connected to the air-line by a } in. pipe about 3 ft. 
long, which was closed by a simple cock. The injection 
of helium was accomplished by opening and closing the 
cock as quickly as possible. The timing drum (see below) 
indicated that helium was flowing into the air-line for a 
period of less than 0-2 sec. in both runs. The quantity 
of helium injected could be calculated from the ambient 
temperature, the volume of the pressure-vessel, the 
initial pressure, and the pressure after the residual gas 
had returned to ambient temperature. The gauge was 
tested after the experiment and found to be accurate, 


Fig. 2. Sampling pipe and branch, with cooling bath. 


Sampling—Samples of the stack-gas were taken through 
a 4 ft. long 1 in. pipe projecting into the stack, This pipe 
was vented freely to the atmosphere, and the pressure- 
drop (9 p.s.i.) was sufficient to cause the gas to flow 
through the pipe at a high velocity. 12 side-branches of 4 in. 
L.D. tube were brazed into the pipe, and passed through 
a water-bath so that their outer ends remained cool 
(Fig. 2). The sample vessels were attached to these 
branch tubes by rubber pressure-tubing. A sample vessel 
is shown in Fig. 3. A three-way tap C gave access either 
to the purge vessel or to the sample-vessel proper. Each 
vessel was connected via a header-pipe to a vacuum 
pump, and before sampling both vessels were evacuated 


Fig. 3. Purge and sample vessels. 


down to tap C, which was turned so as to isolate them 
from the stack-gas and from one another. Tap A was 
then closed (tap D could be left open if the vacuum pump 
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Run 1. 
Nitrogen flow rate : n = 1-85 lb. mole /sec. 
Calculated time of passage from injection-point to base of regenerator = 5-4 sec. 


Time after entry of He t (sec.) | 96 | 151 | 20-6 | 25-2 | 30-7 | 41-8 46-8 518 
(He/N,) 108 108 | | 146 | 351 | 276 | 208 | 160 | 725 55-8 
(He/CO,) « 108 a’ 926 2180 1850 1310 476 374 
e/a’ | - 0-149 0159 0155) 0-151) 0-149 
(He/Ny) x 108 108, a’, 0 | | 116 | 73:5 | 57-7 
L | 


co 
=0155; (Orsat) = 0-148 
alm Ny 


Helium injected ; Q’ = 0-0195 lb. mole ; n | ce dt = 0-0149 Ib. mole 


@o 
N 
Mean residence time : = 83-7 see ; / [eae = 920 
0 0 


Run 2. 


Nitrogen flow-rate n = 1-87 lb. mole /sec. 
Calculated time of passage from injection-point to base of regenerator = 5-3 sec. 


| | | 
Time after entry | | | | 
of He (see.) 47110 (142 162 18-6 212 «23-1 | 317 40-7 60-1 
(He/N,) 10° | 108 | 0 166 670 147 | 198 | 244 256 | 165 465 | 231 
(He/CO,) x 10° 108 a’ O12 435 961 4350 167017301040 323 156 
| | 
COg/Ny | c/a’ i- O154 0152 0-148 0-145 0-146 0-146 
(He/N,) 10° | 108. a’. | 16-6 | 645 148 | 200 | 248 | 256 154 47-9 23-1 


, 


0 
= 0.148 ; (Orsat) = 0-148 
alm Ny 


Helium injected : Q’ = 0-0191 Ib. mole ; n few = 0-0114 Ib. mole. 


@ 
N 
Mean residence time : = 33-5 sec. ; le tlt / few = 33-3 sec. 
n 
0 0 


was kept running). When it was desired to take a sample, | vessel was isolated ; (c) the sampling vessel was opened 
a rapid turn of tap C through some 220° brought about | to the branch-tube and sucked in a sample of stack-gas ; 
the following sequence of events: (a) the purge vessel | (d) the sampling vessel was isolated. Tap B was closed 
was opened to the branch-tube, and sucked through it | at leisure, and the sample removed for analysis. 

several times its own volume of stack-gas ; (b) the purge By this method the branch-tube was flushed out with 
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the gas to be sampled immediately before the sample 
was taken. 


Time Measurements—The time of injection and the times 
at which the samples were taken were recorded by an 
electrically-operated stylus on a rotating drum. The 
stylus drew a horizontal line except when any one of 13 
electrical circuits was closed, when it moved upwards and 
drew a line at a higher level so long as the circuit remained 
closed. The injection of helium caused one circuit to close 
and remain closed so long as the helium flow continued 
(see Fig. 4). The stylus thus recorded the time and 
duration of the injection. As each sample-vessel was 
filled, a thread of mercury rose into the evacuated arm of 
a small capillary-tube manometer (E, Fig. 3). In doing 
so it momentarily made a contact between two platinum 
wires, and the stylus recorded 2 “ blip" on the drum. 
Action at the injection, sampling and timing-points was 
coordinated by field-telephone. 


SPRING FLAP 


Fig. 4. Helium injection timing device. 


Analyses— The samples were analysed with a Metropolitan 
Vickers MS2 Mass Spectrometer. No modifications were 
made to the instrument, which was currently in use for 
routine hydrocarbon gas analysis. However, in order to 
achieve a high sensitivity for the detection of helium 
(0-0005 %) the sample pressure on the high pressure side 
of the introduction leak to the ionization chamber was 
increased to 20 mm.1Hg. 

Some of the samples were contaminated by air which 
leaked into the sample vessels ; a correction was made 
for this by determining the ratio of O, to N, in the samples 
(the stack gas itself contained not more than 0-2 mol % 
Oy, dry basis). 


The accuracy of the reported ratios of helium to nitrogen 
was estimated as : 


true (Ile /N,) 
= + 8) 
reported (He/Ny) 
the error 5 having a maximum value of 0-1 and the 
coefficient « (which was used only in the helium material 
balance) being a constant for the samples in any one 
run. The ratios of the true and reported He /CO, ratios 
were similarly related, « having the same value. It was 
found that at the abnormally high sample pressures, mass 


discriminatory effects made the absolute value of x difficult 
to determine. These discriminatory effects, which were 
absent for heavier molecules (M > 20) at 20 mm.Hg., and 
for helium atoms at lower sample pressures were thought 
to be due either to the onset of hydrodynamic flow, or 
more probably to space-charge effects within the ionisation 
chamber. On the basis of synthetic blends which were 
made up and analysed, « was at first believed to lie between 
0-9 and 1-1. However, in view of the discrepancy in the 
helium balance (see * Interpretation of Results ”) further 
calibrating measurements were made under conditions 
simulating more closely those prevailing during the analysis 
of the experimental samples. The calibrations indicated 
that « had the value of 1-33 for Run 1, and 1-49 for 
Run 2. The former gives a helium balance to about 1%, 
but the discrepancy in the latter is still about 10%. 


INTERPRETATION OF RESULTS 


The tables show the reported ratios of He to 
N,(c’) and He to CO,(a’). The ratio (c’ a’) should 
be equal to the ratio of CO, to N, in the sample, 
and hence should be virtually constant in each 
run and equal to the mean ratio obtained from 
the analyses of the stack-gas carried out before 
and after each run. This provides a useful check 
on the consistency of the analyses. 

Multiplying the He CO, ratios, a’, by the mean 
value (c’/a’),, for the run gives a second series of 
values of the He ‘N, ratios in each run. Since 
these contain independent information, they have 
been given equal weight with the directly-deter- 
mined values of ¢’ in plotting the distribution- 
curves. 

From a knowledge of the dimensions of the 
air-pipe and catalyst-riser, and the volume of the 
catalyst present in the latter, it was possible to 
calculate the delay between the injection of the 
helium and its arrival at the bottom of the 
regenerator. This delay was subtracted from the 
time at which a sample was taken, to obtain the 
value of ¢t, the time between the entry of the 
helium at the bottom of the regenerator and 
taking the sample. The calculation of the delay 
(c. 5sec) was subject to the uncertainties dis- 
cussed below in connection with the calculation 
of the mean residence-time in the regenerator, but 
the error involved is not thought to have been 
great cnough to have any substantial effect on 
the subsequent interpretation of the results. 

If ¢ is the mole ratio of He to N, in the stack-gas 
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at time ¢, and 7 is the total rate of feed of nitrogen 
(which virtually all leaves by the stack), the rate 
of escape of He up the stack at time ¢ is nc. Thus 
the total quantity Q of helium injected is given 
by 


Q=n oa. (2) 


ne 


A fraction — . dt of the helium escapes between 


t and (t + dt) — i.e. this fraction of the helium has 
residence-times between t and (t + dt). Thus the 
mean residence-time of the helium as a whole is 


In order to determine Q and @ from equations 
(2) and (3), e’ and c’t can be plotted against ft, and 
the curves integrated graphically. However, each 
curve has a “ tail’ extending beyond the time 
of the last observation. In order to estimate the 
area under the tail, Inc’ was plotted against t 
for each run, and it was found that the later 
values of c’ in each case fell on a straight line 
with slope 

dine’ 


— = 0.065 sec". 
dt 


Assuming that the tails of the ec’ vs. ¢ curves 
continue to obey this relationship, the relevant 
integrals can be easily computed. (In fact, since 
some helium will find its way into the pores of 
the catalyst, from which it will diffuse out slowly, 
the tail probably becomes flatter than the above 
expression suggests. The above method thus 
probably under-estimates the mean _residence- 
time. However the discrepancy is probably not 
large, as the mean residence-time @ computed 
from equation (3) agrees well with that estimated 
from plant operating data. The good helium 
balance obtained in Run 1 also tends to support 
this view.) 

@ can also be found by dividing the quantity N 


of nitrogen present in the regenerator at any 
instant by the nitrogen feed-rate, n : 
N 

a’ (4) 
N is computed from the void-volume of the 
regenerator and the temperature, pressure and 
composition of the gas. The computation is 
subject to various uncertainties : 

(a) The quantity, and hence the volume, of the 
catalyst in the regenerator cannot be measured 
directly and must be inferred from the pressure 
drop in the catalyst bed. (The pore-volume of 
the catalyst was included in the total void-volume, 
of which it comprised about 5°. 

(b) The pressure fell from about 16 p.s.i.g. at 
the bottom to about 9 p.s.i.g. at the top of the 
regenerator. There was a drop of about 1 p.s.i. 
across the grid. Estimated average values were 
used for the spaces above and below the grid. 

(c) The temperature above the grid was fairly 
uniform (600°C); however, there was a con- 
siderable rise in temperature between the bottom 
of the regenerator and the grid (from about 540 
to 600°C) and an estimated average was used for 
this relatively small volume (5%, of total void 
space). 

(d) No information was available concerning 
the variation of composition of the gas with the 
level in the regenerator (although the work of 
Askins et al [1] suggests that it may be uniform 
above the grid). A uniform value of 71 mole % 
N, (wet basis) throughout the regenerator was 
used in the computation. This corresponds to the 
composition of the stack gas. Since the actual 
concentration of nitrogen could at no point have 
been greater than 79 %, or less than 65%, the 
error introduced by this approximation was 
probably considerably less than 10%. 

In the table of results the value of the mean 
residence-time has been calculated from N /n for 
each run, and can be compared with the value 
obtained by integration according to equation (3). 
The agreement is satisfactory considering the un- 
certainties involved. The values of @ obtained 
from equation (3) are used in subsequent calcula- 
tions. 


On the other hand, there is a considerable 
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discrepancy in each run between Q’, the quantity 
of helium calculated from the pressure and 
volume of the container, and the quantity Q 
calculated from equation (2) and the reported 
helium ratios. In view of the internal consistency 
of each run, the agreement between the residence- 
times calculated by each method, and the agree- 
ment (demonstrated later) between the reduced 
distribution-curves for the two runs, it seems clear 
either that the reported quantity, Q’, of helium 
injected was larger by a factor « (different in each 
run) than the true quantity Q; or that the 
reported helium ratios, c’ and a’, were smaller by 


a factor : (different but constant in each run) 
x 


than the true ratios ¢ and a. For the purpose 
of interpreting the results it makes no difference 
which error has occurred. In either case, 
cQ’ c'Q = «, and equation (2) takes the form 


Q- an fe 
0 


while equation (3) remains valid (putting ec’ for 
ce). a is then found to have the values 1-31 for 
Run 1, 1-67 for Run 2. As explained under 
“ Analyses,” further experiments with the mass 
spectrometer indicated that in fact the reported 


(5) 


(sec) FOR RUN1 


1 
RUN 2 
o 


Fig. 5. C-diagram for Runs 1 and 2. The area under the 


curve between any two ordinates represents the fraction 
of material having residence-times between the corres- 


ponding values of ¢. 


values c’ should be multiplied by a factor of 1.83 
in Run 1, and 1-49 in Run 2. The helium balance 
is thus very good for Run 1; the discrepancy in 
Run 2 remains unexplained. 

It is desirable to plot the results in dimension- 
less form, as a reduced distribution-curve, or 
C-diagram [2]. This leads to a curve which is 
independent of the quantity of tracer used, and 
which directly represents the distribution of 
residence-times. The ordinate, R{t}, is the ratio 
of the concentration of helium in a sample taken 
at time ¢ to that which would be obtained by 
dispersing the helium uniformly throughout the 
regenerator, and was calculated as follows : 


Ne Nac’ Nac’ 
Rit} = ~ => 
an fe dt 


(6) 


— 
| edt | edt 
The abscissa is the ratio of ¢ to the mean residence- 
time : 


(7) 


It will be noticed that no figures derived from 
the plant operating data (e.g., n, N) are involved 
in calculating R and t @ by this method. Only 
the reported values of the He N, (and He/CO,) 
ratios and the measurements of time are used. 
Equations (6) and (7) are exact as long as the 
error « is constant in each run. 

R is plotted against t/@ in Fig. 5 to give the 
C-diagram for each run. Since operating con- 
ditions were almost identical during the two 
runs, the two curves should also be identical In 
fact the points appear to lie on a single curve. 


— fevae 
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This is strong evidence for the internal con- 
sistency of the analyses for each run; it also 
@ 


suggests that the methods of evaluating | ea 
0 
and | e'tdt are not subject to any substantial 
0 


errors. The area under the curve between any 
two values of t @ is equal to the fraction of the 
material which has residence-times between the 
corresponding values of ¢. 

The fraction, F{t}, of the helium (or nitrogen) 
spending less than a specified time, ft, in the 
regenerator, is given by 

t 
= | 0). (8) 


Values of F can therefore be obtained by graphical 
integration of the C-diagram, and are plotted 
against ¢,@ in Fig. 6, which is known as an F- 
diagram (2). 

Discussion 


The most striking feature of the diagrams (Figs. 5 
and 6) is that none of the gas spends less than one 
quarter of the mean residence-time in the regenera- 
tor, and only some 5%, of it spends less than 
half. This proves that under the conditions of 


(sec) FOR RUN 1 
° 8 16 24 32... 40 48 56 


i T T T T 


° ze) 
Yo 
Fig. 6. F-diagram for Runs 1 and 2. The value of F fora 
given value of ¢ is equal to the fraction of material spending 
times less than ¢ in its passage through the system. 


operation the gas in the regenerator (above the 
grid) is not completely mixed; in other words, 
there is no effective back-mixing between the 
top and bottom of the regenerator. The hold-back 


of the system is 0-17, whereas for a completely- 
mixed system it would be 0-37, and for piston-flow 
zero. (The hold-back is a measure of the departure 
from piston-flow, and is equal to the area under 
the F-diagram up to the ordinate = 1; if 
the total quantity N of nitrogen present in the 
regenerator at a given instant were labelled, the 
hold-back would be the fraction of labelled 
nitrogen still present after a quantity N of 
unlabelled nitrogen had flowed in [2]). Only 
about 5°, of the gas spends more than twice the 
mean residence-time in the system. 

It is not possible to deduce from the experi- 
mental results the exact distribution of residence- 
times of the gas in its passage through the 
fluidized bed itself, since little is known a priori 
about the nature of the flow in the space above 
the bed, which has a profound influence on the 
overall distribution as measured. However, a 
surprising amount of information can be extracted 
by the method of analysis to be described. 

It seems probable that the gas in the space 
below the grid is completely mixed, owing to the 
impingement of various streams at relatively high 
velocities. If this is so, the entering helium is 
instantaneously dispersed throughout the inter- 
stitial gas-space below the grid; its mean 
residence-time in this space is calculated (in Run 
1) to have been 2-4sec. The pressure-drop of 
1 p.s.i. across the grid makes it probable that the 
gas flow-rate is uniform through all parts of the 
grid. The high velocity between the bars makes 
it highly unlikely that there is any back-mixing 
across the grid. 

The height of the catalyst bed above the grid 
was estimated, from the pressure gradient and 
total pressure-drop in the bed, to be about 17 ft. 
(This implies a void-fraction of 0-73, excluding 
pores). Hence the total void-volume in the bed 
(including pores) was calculated to be 29% of 
the gas space above the grid, and the mean 
residence-time of the gas in the bed was taken 
to be 29°, of that between grid and stack, or 
8-5 sec. 

If it is further assumed that there is no back- 
mixing between gas which has left the catalyst- 
bed and that within it, the flow through the 
regenerator may be represented as the flow 
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through three compartments which discharge into 
one another in series (Fig. 7). 

Ginuitanpd and Mason [3] showed that the 
residence-time distribution in their model fluid- 
ised beds could be represented approximately by 
expressions of the form : 


(9) 
F 1 =) t 1 ‘ 


where S§ is a parameter which depends mainly on 
the geometry of the bed, and @ the mean resi- 
dence-time in the bed. When S = © the above 
expression reduces to the F-diagram for piston- 
flow, and when S = 1 we have the expression 
for “ complete mixing.” Intermediate values of 
S represent systems with different degrees of 
hold-back between these limits; the hold-back 
in the system is 1,Se. Some F-diagrams with 
different values of S are shown in Fig. 8. 

If it is assumed that the distribution in the 
large fluidised bed in the catalyst regenerator 
also conformed to equation (9), it can be shown 
that the value of S must lie between certain 
limits. Regarding the regenerator as equivalent 
to three compartments in series (Fig. 7) with 
“complete mixing ” in the first (A) and a distri- 
bution corresponding to equation (9) in the 
second (B), it is possible to relate the observed 
overall distribution to the unknown distribution 
in the third compartment (C) by the following 
expression 


F,{t} = F{r} + R{r} 


0, 


[dt (10) 


where +r = '). 

In this expression F,{t} is the fraction of 
material which spends times less than ¢ in its 
passage through the third compartment ; 0), 0, 
@ are the mean residence times in A, B and the 
whole system respectively ; and the number S, 
which determines the distribution in compartment 


B, is unknown, F {7} and R{r} are the values 
of F and R for the system as a whole (taken 
from Figs. 5 and 6 respectively in the present 
case) at a time r = + (S —1)0, 8; the last 
term (in square brackets) is the slope of the C- 
diagram (Fig. 5) for the whole system at time r+. 


=32 sec 


—— 8 Cc 
BOTTOM Cart. above 
CONE BEO 
O, =2 4sec 6, =8'Ssec 0; sec 
Fig. 7. 


Thus, if a value is assigned to S, the corres- 
ponding F,-diagram for the third compartment 
can be constructed from the data for the system 
as a whole. It is found by trial that if § is taken 
to be less than 2, the F-diagram has a region 
of negative slope, which represents a physical 
impossibility. If S is 5 or less the F,-diagram 
displays a point of inflection (i.e. slope has a 
minimum value). This corresponds to a C;- 
diagram with two peaks, and while such behaviour 
is possible for systems in which, for instance, the 
gas flows through two parallel and dissimilar 
channels, it can be confidently assumed that it 
does not occur in the space above the bed. The 
conclusion is that S must be greater than 5 for 
flow through the fluidised bed, the hold-back 
being less than 0-075. In other words (see Fig. 8) 
the flow of gas through the bed does not depart 
very far from piston-flow. 

It should be emphasised that this conclusion, 
in its qualitative form, does not depend on the 
exact validity of the assumptions made, or on 
exact conformity with equation (9); these 
merely serve to simplify the calculations. Pro- 
vided there is no substantial degree of back-mixing 
between the three “ compartments,” and the 
estimate of the void-volume of the bed is approxi- 
mately correct, the experiment can be taken to 
show that the flow through the bed is much closer 
to piston-flow than to complete mixing. 
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CONCLUSIONS 
Although the distribution of residence-times for 
the gas flowing through the fluidised bed cannot 
be deduced exactly from the experiments, it can 
be concluded with some confidence that it is much 
closer to piston-flow than to complete mixing, and 
has a hold-back of less than 0-1, 
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NOTATION 


¢ = true mole ratio of He to Ny in sample 
c’ = reported value of above 
a = true ratio of He to CO, in sample 
a’ = reported value of above 
Q = true quantity of He injected (Ib. mole) 
Q’ = apparent value of above calculated from pressure 
change in container 
a = 
N = nitrogen contained in regenerator (lb. mole) 

n = feed-rate of nitrogen to regenerator (Ib. mole /sec.) 

0 = mean residence time of nitrogen in regenerator (sec.) 

t = time after injection of helium,less calculated time of 
passage from injection-point to bottom of regenera- 
tor (sec.) 

R = Ne/@ (see equation (4)) 
F = fraction of material spending less than specified time 
in passage 
r=t(— (* 


S = parameter in equation (9) 

Subscripts : 

1, 2, 3 refer to gas-spaces below grid, in catalyst bed and 
above catalyst bed respectively. 
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W. G. Bert: Physical Methods in Chemical Analysis, 
Academic Press, Inc. 640 pages. Price $13.50, 1951. 


As in Volume 1, each chapter in this book has been written 
by an expert in the field, and is concerned with a particular 
aspect of the use of physical methods in chemical analysis. 
The first chapter deals with polarography, and gives concise 
accounts of the theory, apparatus, and techniques used. 
Useful tables of depolarisation potentials are included, 
and over forty typical analyses are summarised. The 
chapter ends with a discussion of polarometrice titrations, 
The basic features of polarography are adequately covered, 
and it is therefore unfortunate that only a few references 
to the available literature have been included for those 
who wish to follow the subject further. This is a general 
criticism of many chapters of the book. The next chapter 
covers the theory and technique of conductivity and 
conductometric titration. As a minor criticism one may 
note that the use of ion-exchange resins is not mentioned 
in the section on the preparation of conductivity water. 
The conductometric titration section includes phenols 
and alkaloids as well as inorganic acids and bases and 
titrations which involve precipitations are also reviewed. 
Potentiometric analysis is discussed in an admirable 
manner in the following chapter, and adequate references 
to literature data and reference books are given. 

An excellent account of electrography and _ electro- 
spot testing includes a historical survey, full details, 
with good diagrams, of the equipment, materials, electro- 
lytes, reagents, ete. used, and a good review of uses in 
identification of pure metal surfaces, alloys, and anions. 
Special applications, such as the structure of steel and non- 
ferrous metals, and a section on quantitative applications, 
complete this useful chapter. Adequate accounts of 


magnetic methods of analysis, and of the determination 
of the area of surfaces of solids, are followed by a chapter 
on surface tension measurements. This reviews com- 
petently the capillary rise, maximum bubble pressure, 
drop-weight and pendant drop, and ring methods of 
determination of surface tension, but would have been 
much improved by a detailed section on applications. 

Vacuum techniques of analysis depend first on the 
production and measurement of low pressures, and a good, 
but perhaps over-long, account of this is given. The 
vacuum fusion method for the analysis of gaseous elements 
in metals, and the determination of carbon by low-pressure 
combustion are next considered in detail. The section 
on the micro-analysis of gases at low pressures covers 
only a small part of the field, however, although literature 
references are given to other work. A discussion of 
molecular distillation and a full collection of references 
completes this somewhat unbalanced chapter. A good 
account of analysis by thermal conductivity comes next, 
and is followed by a useful chapter on the measurement 
of radioactivity for tracer applications. This surveys, 
with the aid of excellent diagrams, the photographie and 
electrical methods of detection, the principles of counting, 
and analytical procedures, An equally useful chapter is 
that on statistical analysis, with particular reference to 
spectrography and spectrographic analysis ; numerical 
examples are given. The final chapters surveys briefly 
chromatographic analysis, and is useful as an introduction 
to the subject. 

The book is well printed, is relatively free from major 
errors, and on the whole provides a valuable compendium 
of methods and techniques of physical methods of analysis 
which should be available in every library. 

R. N. 
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